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0^ . Abstract. The correlation evolution of a high luminosity subsaniple 



of the CN0C2 redshift survey is examined. The sample is restricted to 
galaxies for which the k corrected and evolution corrected R luminosity is 
M/j < —20 mag, where M* ~ —20.3 mag. This subsample contains about 

CJ ■ 2300 galaxies. In consort with 13000 galaxies in a similarly defined low 

^j I redshift sample from the Las Campanas Redshift survey we find that the 

comoving correlation can be described as ^{r\z) = (roo/r)'^(l + z)~^^^''> 
with roo = 5.08 ± O.OS/i"^ Mpc, e = 0.02 ± 0.23 and 7 = 1.81 ± 0.03 over 
the z = 0.03 to 0.65 redshift range in a cosmology with 0,m = 0.2, A = 0. 

^ I The measured clustering amplitude, and its evolution, are dependent on 

the adopted cosmology. The evolution rates for rij\,f = 1 and flat low 
density models are e = 0.9 ± 0.3 and e = —0.5 ± 0.2, respectively, with 
roo — 5h~^ Mpc in all cases. 



Introduction 



O ' The measurement of the evolution of galaxy clustering is a direct test of theories 

for the formation of structure and galaxies in the universe. Within CDM-style 
^ ' structure formation theories galaxies form exclusively within dark matter halos 

whose clustering evolution should be slower than the underlying dark matter. 

Clustering evolution is of empirical and pragmatic interest on scales comparable 
K^ , to the size of galaxies themselves, since clustering leads to galaxy-galaxy merging 

Vh ' of gas and stars which is a physical driver of galaxy evolution. On somewhat 

^ larger scales galaxy groups and clusters contain hot gas which exerts evolution 

pressures on galaxies. 
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The theoretical groundwork to interpret the quantitative evolution of dark 
matter clustering and the qualitative trends of galaxy clustering evolution is 
largely in place for hierarchical structure models. Although the details of the 
mass buildup of galaxies and the evolution of their emitted light are far from 
certain at this time, clustering of galaxies depends primarily on the distribution 
of initial density fluctuations on the mass scale of galaxies. N-body simulations 



of ever growing precision and their theoretical analysis (Davis et al. 1985, Ef- 



stathiou et al. 1985 ) have lead to a good semi-analytic understanding of dark 
matter clustering into the nonlinear regime. Normal galaxies, which are known 
to exist near the centers of dark matter halos with velocity dispersions in the 
approximate range of 50 to 250 kms~^, will not have a clustering evolution 



identical to the full dark matter density field (Kaiser 1984). 

The observational measurement of clustering at higher redshifts is still ma- 
turing, with no published survey of the size, scale coverage, or redshift precision 



of the pioneering low redshift CfA survey ( pavis &: Peebles 1983 ). There are 



two redshift surveys extending out to z ~ 1, the Canada France Redshift Survey 
( Lilly et al. 199q ) and the Hawaii K survey ( |Cowie et al. 199q ). Their redshift 



precision is insufficient for kinematic studies, however both quantify the sub' 
stantial evolution in the luminosity function considerable population evolution. 
Measurement of the correlation evolution of the galaxies in these surveys found 
a fairly rapid decline in clustering with redshift ( LeFevre et al. 1996| , parlberg et 



al. 1991] ) although neither analysis took into account the evolution of the lumi- 



nosity function or was able to quantify the effects of the small survey volumes. 



2. The CNOC2 Sample 

The Canadian Network for Observational Cosmology field galaxy redshift survey 
(CN0C2) is designed to investigate nonlinear clustering dynamics and its rela- 
tion to galaxy evolution on scales smaller than approximately 20/i~^ Mpc over 
the 0.1 < z < 0.7 range. There is substantial galaxy evolution over this redshift 
range (Proadhurst, Ellis fc Shanks 198^ , [Ellis et al. 1996| , [Lilly et al. 1995| , [Cowie 



^t al. 199*6 , Lin et al. 1999 ). Investigating the clustering evolution of a population 



requires some basic ability to recover its progenitor hosts at higher redshift. 

The observing and catalogue design procedures are adaptations and exten 
sions of those for the CNOC cluster survey ( Yee, Ellingson fc Carlberg 199£| , Yee 



^t al. 199^ ). The CN0C2 survey is contained in four patches on the sky subtend' 
ing a total sky area of about 1.55 square degrees. The sampled volume is about 
0.5 X 10^/i~'^Mpc , roughly comparable to the low redshift "CfA" survey used 



for similar measurements at low redshift ( Davis fc Peebles 1983| ) which had 1230 



galaxies in the "semi- volume limited" Northern sample from which the correla- 
tion length was derived. The CN0C2 limiting magnitude oi R = 21.5 magnitude 
efficiently covers the redshift range targeted, 0.1 to about 0.7. The survey layout 
for the high luminosity subsample analyzed here is shown in Figure ffl. 

Photometry is obtained in the UBVRI bands, with the R band fixing the 
sample limit at 21.5 mag. The R filter stays redward of the 4000 A break over 
our redshift range. The other bands provide information useful for determining 
appropriate k-corrections and separating galaxies into types of different evolu- 
tionary state ( an issue not considered in this paper). 
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Figure 1. The layout of the high luminosity redshift sample, M^ < 
—20 mag, sample on the sky (left) and in co-moving distance (redshifts 
from 0.1 to 0.65) versus declination in comoving coordinates (right). 



At low redshift we use the Las Campanas Redshift Survey (LCRS) to pro- 
vide a directly compara ble sample. The LCRS is an R band selected survey 
( ^chectman et al. 1996 ) that covers the redshift range 0.033 to 0.15, with R 



band flux between 15.0 and 17.7 mag. The same correlation analysis programs 
were used with for the LCRS and CN0C2 data. The only differences are that 
the LCRS data are not k corrected and is only analyzed for a single cosmological 
model, qq = 0.1. 

2.1. The High Luminosity Subsample 

The evolution of correlations with redshift is sensitive to luminosity evolution. 
That is, usually the goal is to compare the evolution of the correlations of the 
same "primary" galaxies, which have likely evolved in luminosity as a result of 
stellar evolution of the existing stars, star formation, and merging of companion 
galaxies. Over the redshift range discussed here these effects are significant, but 
basically perturbative on a stable underlying population. In the absence of a 
luminosity correction one will be adding in galaxies of lower intrinsic mass with 
increasing redshift. Lower lum inosity galaxies are less correlated ( parlberg et 



|al. 1998 , Loveday et al. 1995 ) which if over-represented in the higher redshift 



will lead to an erroneously rapid correlation evolution. From previous studies we 
concluded that the luminosity function can be approximated as evolving such 
that the absolute magnitude varies as M(z) = M(0) — Qz with Q ~ 1 ( [Lin et 



|al. 1999 ). It should be noted that the this approach is a statistical correction 



designed to identify the same population at two redshifts and will work even 
for an evolving rate of bursting star formation, provided that evolution has 
no significant mass dependence. This assumption should be adequate for our 
high luminosity subsample ( Lilly et al. 1995| , Lin et al. 1999) but may fail for 



at lower luminosities. The volume density of galaxies makes no difference to 
the correlations. For the correlation analysis here we will use galaxies with 
Mj^^ < —20 mag, which defines a volume limited sample over the 0.1 < z < 0.65 
range. The resulting sample (for ^=0.1) contains 2285 galaxies. 

The LCRS data are evolution corrected with the same Q as the CN0C2 
data, although at a mean redshift of about 0.1, this makes very little difference. 
The resulting low redshift sample contains 12467 galaxies that are used in the 
correlation analysis. 



3. Real Space Correlations 

The CN0C2 sample is designed to measure nonlinear clustering on scales of 
lO/i^^ Mpc and less. The clustering is quite naturally measured in terms of the 
two-point correlation function, ^(r), which measures the galaxy density excess 
above a random distribution, no. at distance r from a galaxy, n{r) = no[l-|-^(r)] 
( Peebles 1980| ). Measurement of the real space function ^(r) is not straightfor- 
ward with redshift space data in any case, and yet more difficult with only a few 
hundred galaxies per redshift bin. The projected real space correlation function 
removes the peculiar velocities of redshift space at the cost of making a choice for 
the length of the redshift column over which the summation is done. Noting that 
the correlation function measures the fractional variance in galaxies in different 



volumes we check as much as possible that our evaluation of the correlation does 
not artificially add or remove variance which will bias the result. 

3.1. The Projected Real Space Correlation Function 

The correlation function is a real space quantity, whereas the redshift space sepa- 
ration of two galaxies depends on their peculiar velocities as well as the physical 
separation. Although the peculiar velocities contain much useful information 
about clustering dynamics it is mainly a complication for the study of configu- 
ration space correlations. The peculiar velocities are eliminated by integrating 
over the redshift direction to give the projected correlation function, 



Wp(^p) = / i{\Jrl + rl) dr^ (1) 

J —Rr, 



( Davis &: Peebles 1983 ). li Rp = oo and ^(r) is a power law, ^(r) = {ro/ry 



then Equation | integrates to Wp{rp)/rp = r(l/2)r((7 - l)/2)/r{-f/2){ro/rp^"' 



(Peebles 198C). However, in a practical survey summing over ever the entire 
redshift extent of the survey leads to little increase in the signal and grow- 
ing noise from fluctuations in the field density. Here we are focussed on the 
non-linear correlations, ^ > 1 which suggests we use Rp of several correlation 
lengths. A range of Rp and the resulting errors in rgo (see Equation y below) 
are displayed in Figure ^. We select Rp = Wh~^ Mpc as a conservative choice, 
having the largest (and most stable) errors. Similar results are obtained for 
any 10/i~^ Mpc < Rp < 50/i~^ Mpc, but with increasing statistical fluctuations, 
most notably in the estimated errors of the fltted correlations. 

The choice of a statistical estimator of the correlation function depends on 
the application. With point data the basic procedure is to determine the average 
number of neighboring galaxies within some projected radius and redshift dis- 
t ance Rp. We estimate Wp{rp) using the simplest estimator, Wp = DD/DR — 1 



( Peebles 1980 ), which is accurate for the nonlinear clustering examined here 



and faster than methods which include the RR sum. The DD and DR sums 
extend over all four patches in CN0C2 and the six slices of LCRS, so that patch 
to path variations in the mean volume density become part of the correlation 
signal. The patch to patch variation is used to estimate the error. 

Estimated projected correlation functions, in co-moving co-ordinates using 
Rp = 10h~^ Mpc, are calculated for the LCRS galaxies bounded by redshifts 
[0.033, 0.15] and six redshift bins for the CN0C2 data, [0.10, 0.20, 0.26, 0.35, 
0.40, 0.45, 0.55, 0.65]. The flrst CN0C2 redshift bin is the least populated 
with 185 galaxies and the fourth bin has the most with 602 galaxies. Relatively 
narrow redshift bins helps to reduce any problems associated with the detailed 
shape of the variation of n{z) over the redshift range of the bin. To fit the data 
to a specified function requires error estimates at each point. 

The correlations are fit to the projection of the power law correlation func- 
tion, ^(r) = {ro/r)^ taking the errors to be the 1/yDD, the weighted counts in 
each Tp bin. The fits use only the 0.16 < r^ < 5.0/i~^ Mpc range where there 
are minimal complications from slit crowding and small correlations. All results 
here are derived using co- moving co-ordinates, and normalized to a Hubble con- 
stant Hq = 100/ikms~^ Mpc~^. The results displayed in Figure y are derived 
assuming a background cosmology of Q^m = 0.2, VL\ = 0. 
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Figure 2. The derived roo (see Eq. |2|) as a function of the integration 
length, Rp used to define Wp{rp). The 



confidence intervals are 
shown. The diamonds are plotted when x^ exceeds the formal level 
for 90% confidence. This most likely arises because the variances used 
to calculate x^ are estimated from the dispersion of the four patches, 
which will sometimes lead to erroneously small variances and hence 
large x^ values. 



3.2. Errors of the Correlations 

A straightforward approach to random error estimates is to take advantage of 
our sample being distributed over a number of separate patches. In each patch 
we fit four for CN0C2 and six for LCRS. The open circles are the results for the 
four individual patches CN0C2 fields (the individual LCRS slices are so similar 
that they are not displayed). The solid points give the result from the combined 
data, along with the estimated error. 

As the size of a correlation survey grows there is a systematic increase in 
the derived correlation length. This is is illustrated in Figure 0. The straight 
mean of the CN0C2 rg is 3.2h~^ Mpc, the median is 3.4/i~ Mpc, and the 
mean of Tq is 3.5/i~^ Mpc, whereas the four patches analyzed together find rg 
is nearly 4.0/i~^ Mpc. This raises the question as to whether the correlations 
have converged within the current survey. We note that the expected variation 
from patch to patch for the given volumes with narrow redshift bins is about 
45%, which is consistent the difference between a correlation length of 3.5 and 
4:.3h~^ Mpc. In the combined sample with bigger bins we expect that there 
could be as much about 10% of the variance missing, which would boost the 
correlation lengths another 5%. 
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Figure 3. The correlation lengths (normalized to 7 = 1.8) as a func- 
tion of redshift for Ojvf = 0.2, ^\ = 0. The darkened diamond is from 
the LCRS with g=0.1. The errors are estimated from the variance of 
the four CN0C2 sky patches (shown with plus, asterisk, circle and 
cross symbols for the 0223, 0920, 1417 and 2148 patches, respectively) 
and the six LCRS slices (not shown since they are small). The solid 
points are the result of the combined correlation analysis of the four 
fields and are in general always larger than the mean of the individual 
fields, since they include the field to field variance. 



4. The Evolution of Galaxy Clustering 

The evolution of correlations is conveniently described with the "epsilon" model, 

ro(z)=roo(l + z)-(=^+^-^)/^ (2) 

where all lengths are in co-moving units. The parameter e measures the rate of 
growth of the mean physical density of neighboring galaxies. If e = then there 
is no change in the physical density with redshift. Positive e indicate a decline 
of clustered density with increasing redshift. The derived ro(z) along with the 
estimated errors, a{rQ(z)) are fit to the e model with a formal x^ which reports 
both the suitability of the model and the parameter confidence interval. 

The correlation lengths for the CN0C2 and LCRS analyzed in precisely the 
same way for our standard Rp = 10/i~^ Mpc and Q = 1 are shown in Figure ^ 
It is immediately clear that there is relatively little correlation evolution of this 
sample. It must be borne in mind that the sample is defined to be a similar set of 
galaxies with luminosities larger than about L*, with luminosity evolution cor- 
rected, that accurately approximates a sample of fixed stellar mass with redshift. 
Samples which admit lower luminosity galaxies, or do not correct for evolution, 
or are selected in bluer pass bands where evolutionary effects are larger and less 
certainly corrected, will all tend to have lower correlation amplitudes. 

The e model fits to the measured correlations are shown in Figure Q. The 
results are clearly a strong function of the cosmological world model. However, 
the strength of the correlations is truly impressive for any of the cosmologies. 
In the case of the low density fiat model, the distances to the highest redshift 
galaxies are sufficiently large that the mean luminosity of the highest redshift 
galaxies is larger than those at lower redshifts. The best fit e values are 0.02ib0.23 
for Qm = 0.2, Qa = 0, 

At redshifts beyond 0.1 or so, the choice of cosmological model has a sub- 
stantial effect on the correlation estimates. Relative to high matter density 
cosmological models. Low density and A models tend to have larger distances 
and volumes, which causes the correlations to be enhanced. The LCRS data are 
analyzed only within the qq = 0.1 model. The correlations for three cosmologies, 
flat matter dominated, open, and low density A, are shown in Figure y. The 
evolution rates for the flat and flat low density models are e = 0.9 it 0.3 and 
e = —0.5 lb 0.2, respectively, with vqq ~ 5h~^ Mpc in all cases. These are marked 
with plus signs in Figure ^. 

The effects of alternate values for the luminosity evolution are shown in 
Figure H with crosses indicating the results for Q = and Q = 2, with the 
adopted value being Q = 1. The absolute magnitude limit remains M/j = — 20 
mag in all cases. In the absence of any allowance for luminosity evolution, 
Q = 0, galaxies of lower current epoch luminosity are increasingly included with 
redshift. Since lower luminosity galaxies are generally less correlated that leads 
to an artificially large e, the effect over this redshift range being very roughly 
Ae Pi 0.5AQ. 
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Figure 4. The x^ confidence levels for fits to the e model for the 

^M = 0.2, ^A = model. The contours are for 68%, 90 and 99% 
confidence. The plus signs mark the results for Qm = 1,^A = 0, 
peaking at e ~ 0.6 and ^m = 0.2, Qa = 0.8 peaking at e = —0.4. The 
x's show the results for no evolution correction, Q = 0, in which case 
e = 0.65, and Q = 2, evolution, where e = —0.3. 



5. Conclusions 

The primary conclusion is that correlations evolve very weakly, it at all, with 
redshift. The main correlation results of this paper are contained in Figures ^ 
Out to redshift 0.65 the correlation evolution of high luminosity galaxies, k 
and evolution corrected, can be described with the double power law model, 
^{r\z) = (roo/r)i'(l + z)-^^+'\ with roo = 5.08 ± O.OS/i"! Mpc, e = 0.02 ± 0.23 
and 7 = 1.81 it 0.03. There is no evidence in the current data for a change in 
the slope of the correlation function with redshift. The clustering results are 



consistent with galaxies being located in dark matt er halos ( Colin, Carlberg &: 
Couchman 1997 , Baugh et al. 199q , Colin et al. 1999 ) but are also consistent with 



"light-traces-mass" in low density universes where the mass clustering evolves 
very slowly with redshift. 

Environmental factors clearly play a major role in the development of galax- 
ies, as most clearly seen in the morphology-density relation (pressler 198C, 



Hashimoto et al. 199^ ). We have established here that the neighborhood den- 



sity, at separations between 0.1 and about 5h~^ Mpc, is essentially non-evolving 
with increasing redshift. The implication is that any evolution in the popula- 
tion cannot be attributed to a significant change in the environmental density 
of neighboring galaxies. In a flat low-density model there is a slight rise of the 
mean density of neighbors into the past. The high luminosity luminous galaxies 



examined here do not as a population evolve much over the z < 1 range ( Lilly et| 



al. 1995| , [Cowie et al. 199^ , Lin et al. 1999| ) which could be the outcome of either 



no change in the environment, or, these galaxies being relatively insensitive to 
environmental change. 
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